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Abstract 

An evaluation of the application of the high extra-lunar 
flux in pre-mare times to more general problems of early solar 
system history is attempted by combining the results of dynamic 
studies with lunar chronological data. 

There is a 2 to 4-fold contrast in the integral impact flux 
between the Apollo 14 and 16 sites and the older mare surfaces. 
This is judged insufficient to account for the contrasting li th- 
ology between these two sites: basalts and soil breccias in the 

maria, annealed breccias and impact melts in the highlands. There- 
fore these rocks and thei r 3 . 9-4.0 b.y. ages are thought to predate 
the surfaces in which they are found. Estimation of the flux 
needed to produce these lithologies and difficulties associated 
with extrapolating this further back in lunar history give sup- 
port to the "cataclysm” hypothesis of Tera, Papanastass iou and 
V/asserburg. 

Dynamical studies permit separate evaluation of the possible 
sources for both the "normal" flux during the first 600 m.y. years 
of lunar history as well as the "peak" which apparently occurred 
4.0 b.y. ago. The most likely sources for the normal flux are 
comets from the vicinity of Uranus and Neptune. The most prom- 
ising source for the peak is tidal disruption by Earth or Venus 
of a ceres-size asteroid initially in a Mars crossing orbit. 
Alternative possibilities are suggested. 

Dynamical studies have been carried out in order to deter- 
mine the extent to which a heliocentric flux could be confined 
to the Moon (and Earth). A Monte Carlo method, based on that 
of Arnold (1965) has been used to calculate the relative impact 
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rates of planet-crossing bodies with the rnoon and the terrest- 
rial planets. it Is found that except for nearly circular in- 

itial orbits, the resulting impact density on these bodies is 
similar. Nearly circular initial orbits at the distance of 
Mercury greatly favor impact with that planet. in the case 
of Earth and Venus they tend to share the material with one 
another . 

It is concluded that the time-vari ati on of the flux on these 
planets is closely related to that on .the Moon. 
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1. The extralun^r flux in the early history of the moon. 

In order to understand the evolution of planetary interiors 
and surfaces it is necessary to establish a time scale to which 
the major events in its history may be referred. Prior to ob- 
taining samples of a planet it is possible to establish a rela- 
tive time scale based on geological principles of superposition, 
in which contact relationships between units resulting from 
identifiable events are used to determine which of the events 
occurred first. Unless erosion processes are too severe, it 
is possible to go further and establish a planet-wide time scale 
based on crater frequencies and crater morphology. A time-related 
parameter 



is thereby determined, where $>(t) is the flux of impacting bodies 
at time t in the past. Work in a number of laboratories is es- 
tablishing the relationship between absolute time, t, and the 

87 

flux time, r for the moon by determination of Rb -Sr 07 , K 0 - A r 4 0 
and U,Th->Pb ages of suitable lunar rocks. Under fairly general 
circumstances, as will be explained later, establishment of this 
time scale for the moon represents a major step in establishing 
a cratering time scale for all of the inner planets, as the extra-1 
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flux forms a link connecting the history of the terrestrial planets. 

An early result of the task of relating to t was the dem- 
onstration that was not a linear function of time. Although 
the integrated flux indicated by crater frequencies at the Apollo 11 
site in Mare Tranqui 1 i tat i s was cc.l'that of the lunar highlands 
(Gault, 1970) the age of the basalts flooding this region was found 
to be 3.6xl0 9 yrs, the greater part of the entire 4.6x10 s yr age 
of the moon ( Papanas tass i ou et al , 1970, Turner, 1970). More 
detailed information on the relationship between and t has been 
provided by subsequent lunar missions (Soderblom and Lebofsky, 

1972, Soderblom and Boyce, 1972, Boyce and Dial, 1973). 

Particularly striking is the rapid decline in flux between 
the time of filling of the oldest mare surfaces, and the younger 
pre-mare surfaces (Cayley and Fra Mauro). As will be discussed 
in somewhat greater detail subsequently, there is some uncertainty 
as to the exact time interval over which this decrease took place. 
However it is clear that the decrease in flux during a few hundred 
million years about 4xl0 9 years ago was greater than the decrease 
in flux during all of subsequent lunar history. This has been 
interpreted by Tera et al (1973, 1974) as indicating a lunar 
"ca tael ysm“ , a peak of extralunar bombardment about 3.9xl0 9 years 
ago. 

Both the'Fra Mauro (Apollo 14) and Descartes Highlands 
(Apollo 16) sites are relatively lightly cratered terra surfaces 
which have received 2-4 times the impact flux of the older mare 
surfaces in Mare T ranqu i 1 i tat i s and Mare Serenitatis. At both 
these sites igneous rocks (including probable impact melts) and 
shock metamorphosed breccias have yielded ages of 3.84 to 4.01xl0 9 
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years by the Kb- Sr and K-Ar methods ( Papanas toss i ou and V/asserburg, 
1 97 1 , 1972ab; Comps ton et al , 19/2, Mur thy el al , 197?-, S tot tier 
et al , 1 973, Turner et al , 19/3) and U-Pb measurements indicate 
that this was also a time of U-Pb fractionation, presumably assoc- 
iated with vol at i 1 i zat i on of Pb (Nunes et al , 1973, Tera and 
Wasserburg, 1972ab). Most, if not all, of these ages are best 
interpreted as impact- rel ated. Post-mare impact melted rocks 
and wel 1 -anneal ed breccias are notably absent from the mare sites 
sampled. Although such rocks are undoubtedly produced by the 
same flux which has produced melted rocks on the earth (e.g. 
the Ries basin), the overall effect of this flux for the past 
3.6 b.y. has caused a negligible resetting of ages on the mare 
surfaces. A few post-mare breccias have been investigated from 
the Apollo 15 site (Stettler et al., 1973; Mark et al , 197*0 and 
the age of the mare basalts in these rocks has not been noticeably 
affected by their later breed at Ion. Therefore it is thought 
unlikely tha t i nc reas i ng this integrated flux by a factor of 
only two to four would result in the impact metamorphism dated 
by the Apollo 14 and 16 rocks. Consequently, it is believed that 
these 3 .84-4 .01 xl 0 9 yrs ages predate the surfaces on which they 
have been collected. On. the other hand, Rb-Sr ages of 3-77 to 
3 .82x1 0 9 yr are obtained on Mare Serenltatis basalts at the Apollo 
17 site (Even sen et al , 1973, Tera et al 1974). These mare 
ages nearly overlap those found at the terra sites. This 
smal 1 di fference in ages makes it difficult to calculate pre- 
cisely the rate of decrease of the flux of impacting bodies, 
as the difference between the mare and terra ages are compar- 
able to differences found for the same rock in different 
laboratories by the same method, and also comparable to the 
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uncertainty introduced by half-life errors In comparing ages 
measured by different methods. Probably the most consistent 
way of calculating the change in flux is to use four Rb-Sr ages 
measured at the California Institute of Technology: Fra Mauro 

rocks 14073 and 14310 (3.88+. 03 b.y.), Descartes Highlands rock 
68415 (3.84+. 01 b.y.) and Mare Serenitatis basalt 75055 (3.77+.G6 b.y. 
Use of Rb-Sr ages on the same or related rocks from other labora- 
tories will shift the absolute ages somewhat but will not change 
the differences appreciably. Fitting the measured integral flux 
to these ages, as well as to those of the Apollo 11, 12 and 15 
mare basalts, and the present measured flux requires a component 
of flux decaying with a half-life of 40+10 m.y., joined to at 
least two longer- lived components. 

The gap between the youngest terra rocks and the oldest mare 
basalts becomes even smaller if the 3.9 b.y. mare age found by 
Stettler et al (1974) is included. However it is possible that 
these older mare samples somewhat predate the surface on which 
they are now found. 

The half-life of the short-lived component will be lengthened to 
about 60 m.y. if 14310, 14073, and 68415 are not interpreted as be- 
ing impact melts, but as minor extrusive or intrusive rocks, post- 
dating the surface on which they are found. The meteoritic com- 
ponent found in these rocks (Anders et al, 1973) argues against 
this interpretation. 

The foregoing discussion made no assumptions regarding the 
relationship of the rocks studied to bas i n- formi ng events, e.g. 
the rel a t i onsh i ps of the Apollo 14 samples, the Fra Mauro form- 
ation, and Imbrium ejecta. Such a discussion is not necessary 
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in order to demonstrate the rapid decay of the extra] unar bom- 
bardment. However these results are entirely consistent with the 
view that the Fra Mauro formation is an Inibrium ejecta blanket 
and that the rocks collected near the, rim of Cone Crater are 
derived from the Fra Mauro formation. As discussed above, the 
age of the Fra Mauro surface is known to be less than 3-88 b.y. 
if 14310 and 14073 are impact melts and less than 3-95 ""4 .01 b.y. 
if they are small-scale igneous rocks of internal origin. 

Compston et al (1972) have presented reasons for believing 
that the 3*95 to 4.01xl0 9 year mineral ages found for basaltic 
clasts in 14321 are metamorphic ages, their resetting resulting 
from metamorphism in a hot Imbrium ejecta blanket. Although the 
difference is not far outside of experimental error, this conclu- 
sion is at least nominally inconsistent with the 3.88 b.y. ages 
of 14073 and 14310 being Pre-lmbrium. 

We have carried out some Rb-Sr measurements which show it 
to be unlikely that these ages are metamorphic. We have deter- 
mined a mineral isochron for a basaltic clast in 14321 (Fig. 1) 
as well as for a microbreccia clast in the same rock (Fig. 2). 

The ages calculated from these isochrons are equal, well within 
experimental error. However the initial Sr B7 /Sr GG ratios are 
entirely different, showing that the basaltic clast and micro- 
breccia clast were not in isotopic equi 1 i brium. at that time. 
However this result did not exclude the possibility that the 
basalt did not equilibrate internally, even though such equili- 
bration did not extend between clasts. fn order to check this 
possibility we have obtained Rb-Sr data at the contact between 
a basalt clast and adj acen t microbreccia clasts, all within 1mm 
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of one another (Fig. 3). If the surface region of the basalt 
clast had equilibrated with the adjacent microbreccia clasts, 
all three would lie along an isochron of slope 4.0 b.y. This 
is not the case. The basaltic material does not deviate at all 
from the previously determined internal basalt isochron, and 
there is no suggestion that any of the more radiogenic strontium 
from the microbreccia has entered into the basalt. Therefore 
it seems unlikely that the basalt has reequilibrated completely 
internally, if this rock represents a fragment of, or a clast wi thi n the 
Fra Mauro formation. However, as discussed by Chao (1973) it is 
by no means necessary to believe that the Fra Mauro formation at 
the Apollo 14 site was ever very hot, 1 and Is more likely to rep- 
resent a debris blanket more analogous to the bunte breccia of 
the Ries basin. In this case the 4.0 b.y. ages found for the 
basalt and microbreccia represent Pre-lmbrlum ages, consistent 
with the previous i n te rpretat i on that such ages, as well as those 
slightly younger, predate the cratered surface at this site. 

Following this line of thought, the ubiquitous - 4.0 b.y. 
ages found in the highlands must represent the effect of an in- 
tegrated flux considerably greater than the fairly moderate 2 to 
4-fold increase found at the Apollo 14 and 16 sites. In order 
for these ages to dominate in the manner they do a much greater 
integrated flux, e.g. a further 10-fold increase, is required in 
the interval 3-95 to 4.05 b.y. Alternatively, or possibly equiva- 
lently, this could represent a period of major basin formation. 

In this case it is necessary to assign at least the Imbrium and 
Or i en tale impacts , and possibly that of Crisium and Nectaris as 
well to this short interval of time. If this rapid decline in 
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bombardment by basin and/or crater-formi ng bodies is extrapolated 
back in time an integrated flux more than 1000 times that 
found on the oldest mare surfaces would be expected on a 4.4 b.y. 
surface. Such a rapid increase would greatly reduce the prob- 
ability of still finding 4.2 to 4.3 b.y. ages still preserved, as 
found in several laboratories (Schaeffer and Husain, 1973, Kirsten 
and Horn, 1974* Comps ton and Gray, 1974). In addition, it would 
seen) difficult to confine the radioactive KREEP material to the 
Procel 1 arum- (mbr i um region, as found by the orbiting"^ -spectrometer 
experiments. Rb-Sr data on KREEP-rich rocks show that the orig- 
inal Rb enrichment of this material took place about 4.4 b.y. 
ago, and a continued increase in bombardment beyond that found 
at 4.0 b.y. would be expected to result in more uni form distribution 
of this mater i al . 

In view of the foregoing, the most plausible interpretation 
of the data is that given by Tera et al , that the lunar bombard- 
ment went through a peak about 4.0 b.y. ago, and the rapid decline 
in flux following that time was not merely the "tail" of a contin- 
ually declining flux beginning 4.6 b.y. ago. 

2. Possible sources of the early extra-lunar flux. 

Regardless of whether the actual lunar flux history was the 
late-peaking "ca tael ysmi c M one outlined above, or the more fre- 
quently stated one of continually declining bombardment, some 
source which decays with a short lifetime, i.e. ~ 40-80 m.y. is 
required to be still present 600 m.y. after the formation of the 
solar system. 

A fundamental but incompletely resolved problem is whether or 
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not the bodies constituting this source were in geocentric or 
in heliocentric orbit. The 600 m.y. time delay in the decay 
of the flux is much too long for randomly oriented crossing geo- 
centric orbits. However the lifetime does not rule out special 
geocentric orbits, such as bodies in circular orbits intersected 
by the moon as it moves out to greater geocentric distance as a 
result of tidal friction. There is no special reason to believe 
that clusters of bodies in such circular orbits were likely to 
have ever formed. In fact, current theories of the evolution of 
a geocentric swarm (Kaula and Harris, 1973, Ruskol , I960, 1963, 

1972) Involve formation of ‘‘geocentric pi anetes imal s" much closer 
to the earth. Also if ~ 10 aS g objects were in earth orbit at 
~ 50 earth radii for 600 m.y. it seems surprising that no much 
smaller bodies exist at somewhat greater geocentric distances to- 
day. Never thel es s, unders tand I ng of the formative stages of the 
solar system is Insufficient to preclude this possibility. 

The principal difference between impact by bodies In geo- 
centric or heliocentric orbit will be their impact velocity: 

~ 3 Km/sec for geocentric bodies and S 1$ Km/sec for heliocentric 
bodies of sufficient lifetime. The kinetic energy of the impacting 
geocentric bodies is barely sufficient to melt the Impacting body 
itself. At the velocities stated above the kinetic energy per 
gram will be about 25 times as great for the heliocentric bodies 
and much more extensive shock-melting would be anticipated. 

The arguments presented in the preceding section implicitly 
assumed that the energy of the impacting bodies in early lunar 
history was the same as that in more recent times (at present 
30 Km/sec). If the earlier bodies were less energetic, the 



increase in flux necessary to explain the widespread shock meta- 
morphism 4 b.y. ago would be much greater. However in view of 
the special orbital circumstances already required for the geo- 
centric source, it does not seem too much of an additional burden 
to require the bodies to be 100 times larger. 

An even lower-energy source of lunar impacts, moon-orbiting 
bodies ~ 10 23 g In mass has been proposed by Reid (1973)* 

It seems difficult to believe that there i s no way to di sting 
uish between the highland breccias and craters being produced 
by very large slow bodies on the one hand or by smaller bodies 
with — 100 times as much energy per unit mass on the other. How- 
ever this seems to be the case at present, and it would represent 
a major contribution to lunar science if this uncertainty could be 
el imi nated. 

With regard to bodies in heliocentric orbits, there are a 
number of classes of bodies in various kinds of orbits which 
have undoubtedly played a role in bombarding the moon. It would 
be a mistake to seek to explain all lunar impact phenomena in 
terms of a single source, even though 'such interpretations may 
appear more simple. 

In discussing the contributions of various sources it has 
been conventional to speak of various "components 11 of the flux 
which begin to decay 4.6x10° years ago with various half-lives. 
This approach Is all right provided It is not taken too seriously. 
Although an exponential decay is a fair approximation for some 
classes of initial orbits, there are others which provide a flux 
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with a more complex time history. An example of this are Mars- 
crossing bodies with perihelion near Mars. For such bodies a 
discrete interval of. 5 to 1.0 b.y. is required for any members 
of an initial ensemble or "swarm" of such bodies to become 
earth-cross i ng . Therefore the rate of earth and moon impacts 
from such a source will remain zero for a long time, increase 
to a maximum, and then decay with a half-life of ~ 1 0 9 yrs. 

The orbital evolution of a number of classes of initial 
heliocentric orbits have been studied by Monte Carlo iteration 
of Opik's (1951) collision theory ( A rnol d, 1 965 ; Anders and Arnold, 
1965 ; Wetherill and Williams, 1968; Wetherill, 1971 > 1972). 

These results will be used in the subsequent discussion. Although 
there are special cases where this theory i s not applicable, 
particularly where commensurab \ 1 1 ti es in period occur, the gen- 
eral lifetime and history of most planet-crossing initial or- 
bits are treated sufficiently well by this method. In any case 
these calculations are to be preferred over those in which 6pik's 
collision theory is used'without iteration. This is particularly 
true when more than one planet’s orbit is crossed. 

Except for the problem of the probable "cataclysm” at 4 b.y. 
to be discussed later, there is no special difficulty in des- 
cribing plausible sources for the generally declining flux over 
most of lunar history. During the earliest accretionary phase 
of the moon low velocity heliocentric (or geocentric) orbits 
arc required. After the moon grew to near its present size, the 
bulk of any residual ea rth-c ross i ng bodies will be swept up on a 
time scale of ~ 1 0 7 years. A similar time scale will apply to 
bodies with aphelion near Jupiter in orbits similar to those pro- 
viding the source of most of the large meteoroid flux at the 
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earth today. Bombardment by these bodies probably played a 
major role in establishing the present size distribution in 
the asteroid belt and may have been important in supplying 
material for the final stages of the earth and moon. However 
it is doubtful if any record remains of the impact of these 
transient sources. 

On the more relevant time scale of ~ 10 a years there are 
several sources which need be considered. Some members of 
the short-lived initial earth-crossing population will be per- 
turbed into orbits of high inclination, and thei r 1 i fet ime will 
be lengthened to ~2xlO e years. On this same time scale Mars 
perturbations will transfer into earth-c ross i ng orbit Mars 
crossers in initial orbits similar to Eros, i.e. those with 
initial perihelion near Earth. At the present time both of 
these sources are of such minor importance that it is very dif- 
ficult to estimate their initial strength. The total flux re- 
quired to impact the moon during the first ~ 500 m.y. after the 
end of terminal accretion may be estimated at ~10 22 -10 23 g. Even 
for bodies with a high probability of terminating their helio- 
centric history by planetary impact, rather than by ejection 
from the solar system, the efficiency of the moon for capturing 
such bodies is small because of its small cross-section relative 
to Earth and Venus. An adequate source will therefore require 
an initial mass of ~10 24 to 10 3G g mass in the proper orbits. This is 
similar toor greater than the present mass of the entire asteroid 
belt. There is no evidence that such a mass of material was not 
present in these orbits in the early solar system. With regard 
to the "Eros- type*' orbits a somewhat special distribution of 
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initial orbits is required to avoid a residue of very long-lived 
(~2x10 9 yr) Mars crossers greater than the number actually observed 
at present. 

As mentioned before, data from the Pra i r i e Network (McCrosky, 
1967) shows that most of the present flux of large meteoroids 
at the earth is of cometary origin (Wetherill, 197 * 0 , primarily 
non-volatile residues of short period comets with aphelion near 
Jupiter. Therefore it is appropriate to inquire as to the role 
of this source in the early solar system. At the present time 
it is believed that comets are perturbed by passing stars into 
the inner solar system from the Oort cloud of comets at distances 
of 1 0 4 to 10 6 A . U . It is usually thought that the Oort cloud 
was initially generated in the region of the major planets and 
perturbed into its present position by perturbations of the 
major planets and later by nearby stars. This process of populat- 
ing the Oort cloud is roughly the inverse of the process by 
which they are captured into "Jupiter's family" at the present 
time. Those cometary bodies which were originally in the vic- 
inity of Jupiter will be lost from the solar system proper on 
a very short time scale (^1 0 s yr) and will contribute only to 
the terminal bombardment. However those comets formed in the 
vicinity of Neptune and Uranus will evolve more slowly. First, 
perturbations by Neptune and/or Uranus will cause their perihelion 
to cross Saturn. In turn Saturn perturbations will lead to 
Jupiter crossing. Once in Jupiter crossing, the further evolution 
will be similar to those of a modern short period comet, and a 
similar f rac t i on will be b r ought into ea rth-crossing orbits. 

The orbital evolution of bodies initially in Neptune-c \ oss 
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i ng or Uranus crossing has been calculated by the Monte Carlo 
technique of Arnold. The initial orbits chosen are: 


Neptune crossing: 


Uranus crossing: 


a e i 

36.0 A.U. 0.2 6° 

28.0 A.U. 0.15 6° 

30.2 A.U. 0.05 1 t° 


23-5 A.U. 0.2 6° 

18.0 A.U. 0.2 6° 

19.0 A.U. 0.05 4 • 5* 


These calculations were made using the original Arnold pro- 
cedure, rather than that I have used more recently, wherein the 
effect of the free oscillations of the secular perturbations was 
included. At present there are no adequate data on this phenom- 
enon for orbits in the outer solar system. 

Histograms showing the time at which the body first becomes 
Jupi ter-cross i ng are presented in Figs. 4-6. Following Jupiter- 
crossing the orbital evolution will follow the course previously 
calculated for short-period comets. The time scale for this 
subsequent evolution will be ~10 s years, usually culminating in 
ejection by Jupiter into hyperbolic hel iocentric orbi ts . If 
earth or moon impact occurs, i t wi 1 1 be within ?0 7 years of 
initial crossing. Therefore on the time scale of lO^-lO 0 years, 
the time distribution of Jupi ter-cross i ng is equivalent to the 
time distribution of earth or moon impacts. 

The shaded portion of the histograms shown in Fig. 4 and 5 
represent those initial orbits with low eccentricity. The gen- 
eral similarity of the low and moderate eccentricity results 
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shows that the results are not strongly dependent on details of 
the ini t i al orbi ts . 

The cumulative impact flux from these sources is shown (on 
a logarithmic scale) in Fig. 7. There is considerable similarity 
between this curve and the lunar impact rate inferred by various 
authors on the basis of crater investigations (e.g. Boyce and 

Dial, 1973), 

No increase is found at 4 b.y., corresponding to the cataclysm 
of Tera and Wasserburg. The flux of these bodies on the earth 
may be estimated to be: 

present cometary flux _ lifetjme near Neptune x No. ^ in Oort cl oud, 
primordial cometary fTux lifetime in Oort cloud No, in Neptune Region 

The lifetimes are ~10 s and ^,1 Q 1 ° years respectively, and the ratio 
of the number of comets in the Oort cloud to the number in the 
vicinity of Neptune is the efficiency for population of the cloud 
which may be estimated to be -1%. Therefore it is estimated that 
the early cometary flux was ~10 3 4 times the present flux. This is 
adequate to provide the integrated early cratering of the lunar 
terra, and it is quite possible that cometary bodies and their 
inactive cores have been the primary source of lunar bombardment 
over almost all of solar system history. 

3- The lunar cataclysm. 

-Plausible natural sources have been found for the principal 
features of the t irne-dependen t lunar meteoroid flux, with the 
exception of the peak in this flux which appears to have occurred 
4 b.y. ago. 

Some suggestions can be made as to how this missing peak can 
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be supplied from a heliocentric source in as natural a way as 
poss i bl e. 

The 600 m.y. interval between the origin of the solar system 
and the appearance of this peak suggests that it be associated 
with the disruption of a large body belonging to one of the pop- 
ulations previously identified as having a characteristic life- 
time of hundreds of millions of years. The rapid decay of the 
disrupted material must occur from orbits with characteristically 
short lifetimes of <10 e years, such as earth-crossing or with 
aphelion near Jupiter. 

The mass of material impacting the moon in this event must 
be — 1 0 s s g - As in the case of the integrated "normal 11 flux, the 
small cross-section of the moon requires that the total mass dis- 
rupted in the inner solar system be much greater, ~!0 34 g for bodies with 
a high (~50%) probabi lity of impacting a terrestrial planet and 
even larger for bodies with a high probability of being deflected 
into Jupi ter-cross i ng and from thence into hyperbolic ejection 
orbi ts . 

Three heliocentric sources of the necessary 1 i fet ime have 
been i dent i f i ed : 

(1) High inclination (e.g. 45°) Earth-crossing orbit 

(2) Mars-cross ? ng with per i he! i on well within the orbit 
of Mars (e.g. Eros) 

(3) "Comets" in Uranus or Jupiter orbits. 

The latter two sources have the characteristic that although 
Lhe lifetime of the source region is maintained in the ~10 R yr range 
by its initial conditions, d i srupt i on can occur in low inclination 
earth-crossing orbit, leading to a short residual lifetime for 

the Fragments. 
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The only portion of the solar system definitely known to 
contain objects anywhere near the proper mass is the asteroid 
belt. However the required mass of ~10 3 -g is about equal to the 
mass of the largest asteroid, Ceres. There is no special reason 
why or why not a similar body should have been in a suitable 
Mars-cross i ng orbit after the formation of the planets. On a 
time scale appropriate to the observed 600 m.y. interval this 
body will be deflected by Mars into Earth-cross i ng. Earth per- 
turbations will usually cause it to become Venus -cross \ ng as 
well within ~10 6 years. It is necessary that this body have a 
high probability of being disrupted while in Earth-crossing orbit 
rather than impacting a planet prior to disruption. The reason 
for this is that in the latter case it would have only a small 
probability of hitting the moon. In order to avoid an ad hoc 
assumption that the one such body surviving until 4 b.y. took 
the improbable course of hitting the moon rather than Earth or 
Venus, it would be necessary to assume there were many such 
bodies. This would then lead to the problem that the record 
of their subsequent impact cannot be found on the earth. In 
order to produce more than one mare basin the even more difficult 
problem of requiring more than one moon- i mpac t i ng body would 
be incurred. A break-up avoids these difficulties by spreading 
the product of a single impact among the various terrestrial 
planets and the moon at a single time. 

The most obvious way to cause a disruption is for the 1 0 34 g 
body to impact a smaller asteroid (e.g. 10 2o g) wh i 1 e traversing 
the asteroid belt. The difficulty with this is that again an 
improbable event will have been invoked. Calculations of the 
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lifetimes of asteroids (e.g. Wetherill, 1967) show that the 
lifetimes of asteroids of this size are -'lO 11 years, and the 
probability of being destroyed by collision during the ~10 7 years 
it is in Earth-crossing is very small (^ICf*). The somewhat great- 
er density of bodies in the asteroid belt 4 b.y. ago may increase 
this probability slightly. However by this time most of the 
short-lived bodies will have been removed and the effect will 
be sma! 1 . 

There exists a much more probable cause of disruption. The 
Earth or Venus impacting probability of an Earth or Venus cross- 
ing body in an orbit derived from Mars crossing is about 30%. 

For small bodies, about this fraction will end their course by 
hitting the Earth or Venus; the remainder will strike Mars or 
Mercury, make a very close approach to the sun, or be perturbed 
by Earth or Venus into Jupiter crossing and subsequent ejection. 
However prior to impacting Earth or Venus there is a high prob- 
ability that the body will come within ~ 3 planetary radii of 
one of these planets, i .e. within the Roche limit. This i s more 
probable than actually Impacting the planet in the ratio /R^ 3 

where R n is the Roche distance and R is the planetary radius. 

R p r . 

Therefore for bodies of sufficiently low strength relative to 
the force of tidal disruption a much more probable event is 
tidal fragmentation. Most of the remaining fragments will con- 
tinue in heliocentric orbits. Thei r si i ght 1 y different veloc- 
ities following disruption will result in different subsequent 
planetary approaches and to divergent orbital evolution. A 
portion of the fragments will hit the moon, in accordance with 
the relative impact probabilities described in the next section. 

Tliis is a probable course of events leading to a lunar cat a- 
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clysm provided a sufficiently large Mars-crossi ng body were in- 
i t i al 1 y p resent . 

An identical disruption history could be described for an 
unusually large comet derived from Neptune or Uranus crossing. 

)n contrast to the asteroidal situation no comets of sufficient 
size are known. The cometary contribution will be dominated 
by those few bodies which chance to evolve into Earth-crossing 
orbits with aphelion St. 2 A.U. Such bodies will have about a 
5% probability of impacting Earth or Venus, and about a .15% 
probability of hitting the moon. Therefore a comet with the 
apparently improbable mass of ~10 ss g is required. This may be 
compared with the ~10 18 g mass usually assumed for comets. How- 
ever comets have been observed for only a few thousand years, 
and a comet this large would not necessarily be a conspicuous 
object if its perihelion were beyond 2 A.U. There is at present 
one body of ~10 27 g mass in a Neptune-cross i ng orbit (Pluto) and 

the possibility of there being many 1 0 s n g objects there early 

in the history of the solar system cannot be ruled out a_ priori . 

The high inclination of Uranus' rotation axis may be a further 
indication of the former presence of large bodies in this region. 

A conceivable way to produce the cataclysm without invoking 
an unusually large comet would be to assume that the formation of 
Uranus and Neptune were delayed by 600 m.y.. Such a delay, owing to 
their slow rate of accretion, is in itself plausible, as discussed by 
Levin (1972). in this case the cataclysm would correspond to the in- 
itial "spike" produced from Uranus crossers during the first 10 e years. 
This alternative is not very attractive. The effect is not large 
enough. furthermore if Uranus were delayed 600 m.y. It is nat- 
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ural to suppose Neptune would be delayed much longer. This 
would lead to a second “Neptune cataclysm 11 late In solar 
system history, which is not observed on either Earth or Moon. 

Undoubtedly there are other possibilities, but those I 
have thought of present more difficulties than those mentioned 
above. For example the mechanisms involving as teroi d- Jup i ter 
resonances proposed by Williams (1973) and Zimmerman and Wetherill 
(1973) a 1 most ce r ta i nl y contribute to the lunar meteoroid flux. 
However in order to provide the peak in flux 4 b.y. ago a rw10 35 g 
asteroid would have to be totally disrupted In the asteroid belt. 

As discussed earlier, this is an unlikely event. 

Of those processes considered, the most plausible to produce 
this peak is probably tidal disruption of a large asteroid per- 
turbed by Mars into an earth-crossing orbit. The similar pro- 
cess with a cometary source, as well as those involving geocentric 
bodies may be considered as possible al ternat ? ves . 

4. Cratering on the terrestrial planets. 

Now that a reasonable understanding exists concerning the 
absolute cratering time scale of the moon, it is of interest to 
learn what this can tell us about the cratering time scale of 
the other terrestrial planets. This possibility is largely con- 
fined to the heliocentric component of the cratering flux, again 
emphasizing the importance of finding some way to separate the 
effect of these two sources. 

An object in earth-(and moon) crossing orbit will make 
close approaches to the earth much more often than it strikes 
the earth. These close approaches will change the orbital elements 
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of the body. As a result of these changes the body can, for 
example, become Venus crossing, after which close, approaches 
and impacts with Venus become possible also. Continuation of 
this process can lead to the possibility of the body striking 
any of the terrestrial planets. The relative probability of 
striking the various planets will depend upon the initial orbit 
of the body . 

The orbital evolution of various initial orbits has been 
calculated using the development of Arnolds' Monte Carlo method 
described previously (Wetherill, 1969). Some modifications in 
the calculations were made in order that a statistically signif- 
icant number of impacts with the smaller bodies, i.e. Moon and 
Mercury could be obtained. This was done in two ways: 

(1) As described by Arnold (1965), the close approaches 
are considered as random impacts on a target circle centered on 
the planet and oriented perpend i cul a r to the velocity of the pass- 
ing body. The radius of the target circle is chosen to be a 
specified number of planetary radii, K. Impacts on the outer 
half of this target circle, i.e. approaches at distances greater 
than K/2 planetary radii, are weighted differently than those 
on the inner half, in order to include the effect of more distant 
approaches within the sphere of influence of the perturbing 
planet, but beyond the radius of the target circle. This causes 
the perturbing effect of a particular planet to be independent 
of the radius chosen for the target circle associated with that 
planet, provided it is chosen to be large enough to eliminate 
large statistical fluctuations. In order to increase the number 
of interactions with the smaller planets (Moon, Mars, and Mercury) 
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a larger target radius, K, is chosen for these planets. This in 
itself does not increase the probability of impacting the planet 
at all, nor does it increase the perturbing effect of the planet. 

It merely smoothes out statistical fluctuations which would 
otherwise occur if a small planet had to compete on an even 
basis with large planets in the choice of the next planet to be ap 
proached. 

(2) For intersections on the inner half of the target circle 
the only difference between an impact and a close approach is 
that for an impact a random number is found to have a value be- 
tween 0 and ^/K 2 > whereas a close approach occurs when the number 
i s between '♦/K 2 and 1.0. For example, for K =10 planetary radii, 
the area of the inner half of the circle is 25 square planetary 
radii. The probability of impact is V?,b, which corresponds to 
the random number being between 0 and .04. However this choice 
of the range of the random number designated as "impact" rather 
than "close approach" is arbitrary. Therefore close approaches 
resulting from the number being between .04 and .08 would have 
been impacts if impacts had been assigned to random numbers in 
this interval. In this way close approaches, in the inner half 
of the target circle may be "scored" as impacts, and the evolu- 
tion calculation permitted to continue. In order to do this 
properly it is necessary to permit each range of .04 to be 
counted as an impact only once, and also to continue the orbital 
evolution when the random number is between 0 and .04. When 
this procedure of scoring close approaches as impacts is com- 
bined with the procedure described in (1) of augmenting K for 
the smaller planets, the apparent relative impact rate on the 
smaller planets increases, thereby providing a more statistically 
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significant number of impacts. Effectively, the number of runs 
has been increased by a factor (K/2) 3 as each run is now equiv- 
alent to (K/2) 2 runs with different choices of the random number 
r eg i on designated as impact for each run. At the end, the effect 
of obtaining more impacts on the smaller planets is corrected 
for exactly by dividing the number of impacts on each planet by 
the value of K 3 used for that planet. Without these modifications, 
the cost of obtaining statistically significant data for lunar 
impacts would be prohibitive. 

The results of these calculations are shown in table 1. AH 
these probabilities are relative to the Moon, which is assigned 
the val ue of unity. 

The first group (A) of initial orbits correspond to initial 
orbits which penetrate deeply into the inner solar system. All 
have aphelia in the asteroid belt. The first two have perihelia 
within the orbit of Mercury, while the third has perihelion 
within the orbit of Venus. It is seen that the probability per 
unit area of impacting Mercury, Venus, Earth, or Moon are com- 
parable, whereas that of h i tt i ng Mars is somewhat less. The Mercury- 
crossing bodies have a higher probability per unit area of 
striking Mercury than the other planets. The statistical un- 
certainty in these numbers is about 10%. 

The first entry in Set B corresponds to initial orbits sim- 
ilar to those from which most meteoroids become earth-crossing 
at present (Wetherill, 1971)* Initially the earth is at Its 
perihelion, and the geocentric velocity low. This causes the 
gravitational radius of the earth to be large, increasing its 
impact rate relative to the moon. Further orbital evolution 
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results in frequent Venus impacts and an appreciable rate of 
Mercury impact. To some extent Mercury is shielded by Earth 
and Venus and the rate of Mercury impact is about an order of 
magnitude lower than found for the Mercury-crossers in Set A. 

The rate of impact on the Moon and Mars are comparable, as dis- 
cussed elsewhere (Wetherill, 1974). The second entry in Set B repre- 
sents an initial orbit beyond the orbit of Earth and with aphelion 
near Jupiter. Jupiter pertu rbat 1 ons will frequently bring such 
bodies into Earth crossing. Again, the impact rate per unit area 
on all the terrestrial planets is the same order of magnitude. 

A more extreme situation is illustrated by the orbits of 
Set C. Here the initial orbits crossed only one planet. If the 
planets were in perfectly circular orbits no other planet would 
be crossed, because of the conservation of relative velocity in 
subsequent encounters. This is equivalent to the constancy of 
the Jacobi integral in the restricted 3-body problem. In spite 
of this tendency for the orbits to remain in the vicinity of one 
planet, they nevertheless evolve in such a way so as to have a 
comparable probability of impacting all the terrestrial planets, 
again per unit area. A marked exception is the Mercury-crossing 
case, which never becomes Venus -c rossi ng . The small mass of 
Mercury and the large energy change necessary to become Venus- 
crosslng preclude this occurring at a significant rate. 

These results do not include the effect of the mean energy 
of the impacts varying from planet to planet. The energy of 
each impact was also computed. This effect will be particularly 
important for Mercury, because typical impact velocities on that 
planet are about twice those on the moon. 


Therefore a smaller 
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mass will produce a crater of the same size. The results of a 
typical calculation of the velocity distribution are shown in 
figure 8. Depending on the exact size distribution of the impact- 
ing bodies, this will increase the integral number of craters 
greater than a given diameter by a factor of 2 to 3. 

It can be concluded that the heliocentric flux producing 
lunar craters will also produce a similar number of craters per 
unit area on all the terrestrial planets. If such a flux causes 

a “cataclysm" at 3.9 b.y. on the moon, it will cause a cataclysm 

on all the planets. Therefore a reasonable first approxirnat i on 
might be to assume that, except for a constant factor $(t) is 
the same for all the terrestrial planets. There is still a ser- 
ious uncertainty in the constant factor. At the present time 

the meteoroid flux in the lOOg-IO 6 g range as given by McCrosky 
(1967) is dominated by orbits similar to the earth-grazi ng orbit 
Al and the Venus -cross i ng orbit B3. Including the velocity effect 
at Mercury, it is probably most reasonable to use the same pro- 
portionality factor for Mercury and the Moon. Mars is a somewhat 
different case. In addition to the sources prevailing at Earth's 
orbit, Mars will receive impacts from asteroids on the inner edge 
of the asteroid belt. In addition, there is significant transport 
of dust and associated erosion on Mars. Considering these prob- 
lems, as discussed elsewhere (Wetherill, 1974), it Is probably 
best to use a Mars flux/Moon flux factor of 1 to 2 rather than 
the value of 10 frequently used. 

There is still much to be done before an entirely satisfactory 
time scale exists for the terrestrial planets other than the earth. 
However it may be expected that present lunar and planetary studies 
will be a very significant aid in this endeavor. 
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Table 1. Relative Impact Probability per Unit Area 

Moon = 1.0 

A. At least Mars, Earth, or Venus Crossing 


Obj ect Peri hel i on 

Aphel I on 

Earth 
Impact 
Eff . (%) 

Merc. 

Venus 

Earth 

Mars 

P . Encke 

0.34 

4.18 

6.0 

4.21 

2.37 

1 .23 

0.37 

Icarus 

0.19 

1.97 

18.5 

7.50 

1 .97 

1 .08 

0.45 

1959LM 

0.70 

3.61 

5-2 

1 .60 

2.34 

1.34 

r> 

CM 

o 

B . Aphel 1 

on near 4. 

A.U. Earth or Mars Crossing 




1 .01 

4.00 

24.0 

0.30 

0. 94 

2.26 

1.15 


1 .36 

4.67 

0.11 

0.69 

1.26 

1 .78 

0.71 

C. Nearly 

Ci rcul a r 

Initial Orbi ts 





Near: 








Mars 

1 .27 

1 .89 

27.8 

1.08 

2.58 

2.82 

11.9 

Earth 

0.98 

0.99 

45.0 

1 .05 

3.02 

3.99 

0.96 

Venus 

0.67 

0.82 

29.6 

1 .08 

2.44 

1 .36 

0.18 

Mercu ry 

0.33 

0.41 

<.01 

>1.5 

xlO 4 

, - 

- 

- 
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Figure Captions 

Mineral isochron for basaltic clast from breccia 
14321 

Mineral isochron for microbreccia clast from breccia 
14321. The line marked 4.05+ .06 b.y. is the isochron 
drawn through density fractions of clast 184-53. The 

other are "whol e- rock" microbreccia clasts. 

\ 

, Rb-Sr evolution diagram showing lack of isotopic • 
equilibrium at the contact between a basalt ciast 
and adjacent microbreccia clasts 

Histogram showing the time distribution of lunar im- 
pacts resulting from comets initially in Neptune- 
c ros sing orb i ts 

Histogram showing the time distribution of lunar im- 
pacts resulting from comets initially in Uranus-crossing 
orbi ts 

The time distribution of lunar impacts from the 
combined effect of Neptune and Uranus-crossing comets 

Cumulative lunar impact flux from Neptune and Uranus 


sou rces 
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Figure 8. Calculated Impact velocities on Mercury and Moon for 
a body initially in an orbit similar to that from 
which me teor i tes are being derived at present. Al- 
though the absolute values of the velocities are 
somewhat dependent on the initial orbit, the ratio 
of velocities remain within a factor of ^2 for most 


r n i 1 1 a 1 orb i ts . 
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